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FIG. 1. The ordinary decay of a b quark, via b —* cW . 
INTRODUCTION 

While hadrons containing bottom quarks decay weakly, and hence are quite long-lived, all 
beautiful things must one day come to an end. For the typical B meson, the end comes after 
about 1.5 ps. While this provides enough time for the meson to pass through a measurable 
distance within a detector, given today's silicon technology, their lifetime is still so short 
that B mesons can be studied experimentally only by the careful examination of their decay 
products. Hence the study of the bottom quark is essentially the study of its decays. 

It is believed that almost all bottom quarks decay weakly into charm quarks, via the 
PF-emission process depicted schematically in Fig. |l|. Rare b decays, then, are those which 
do not include the release of a c quark into the final state. These may include both Cabibbo- 
suppressed decays, such as those mediated by the transition b — > uW~, and flavor-changing 
neutral decays, such as penguin- induced transitions. While the dominant decay mode of the 
b quark is believed to be well-understood, it is hoped that the rare decays may provide a 
window onto new physics beyond the standard model. Not only may one test the standard 
model by comparing the small predicted rates for rare channels to experiment, but the 
very fact that the charmless channels are suppressed makes them ideal places to look for 
anomalous enhancements coming from new particles and interactions at high energy scales. 

The theory of rare b decays has two distinct parts, which are separated from each other 
conceptually and practically by their dependence on physics at very different energy scales. 
From the "high-energy" viewpoint, rare b decays are mediated by intermediate particles of 
large virtuality, and the challenge is to understand the structure of the quark-level transitions 
which such virtual particles can induce. From the "low-energy" viewpoint, rare b decays 
are mediated by local and nonrenormalizable point interactions, with coefficients which 
are determined at high energies, but at low energies may be viewed simply as coupling 
constants of the theory. The relation between the high-energy and low-energy viewpoints is 
demonstrated schematically in Fig. |^ for two typical transitions. 

From the low-energy viewpoint, the theoretical challenge is to relate the strengths of 
the suppressed nonrenormalizable quark-level couplings to physical properties of observable 
hadrons such as B and The situation is complicated by the long-distance effects of 
the strong QCD interactions. Typically, the structure of bottom hadrons cannot be com- 
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FIG. 2. Rare decays of a b quark, from the high-energy and low-energy viewpoints. 

puted from first principles, and one must find techniques which minimize one's sensitivity 
to uncomputable low-energy effects, while allowing one to extract from experiment as much 
information as possible about high-energy physics. At low energies, what one would like to 
measure experimentally are the coefficients of the nonrenormalizable operators such as those 
pictured in Fig. [| It is the goal of low-energy high-energy physics to make this possible. 

In what follows, I shall review the theory of rare b decays both from the high-energy and 
the low-energy points of view. In contrast to the spirit of the rest of this conference, however, 
my emphasis will be on the physics at low energies. In focusing on these possibly less-familiar 
effects, I hope to convince this "high energy" audience of the important limitations which 
low energy strong interactions place on understanding the physical manifestations of virtual 
high energy interactions. The good news is that much work is still in progress to minimize 
these limitations and to maximize the fundamental discovery potential of experimental b 
physics. 

HIGH ENERGY VIEWPOINT 

The decays of b quarks, both ordinary and rare, are generated by virtual interactions at 
some high scale M 3> m^. At lower scales /i < M, these interactions generate nonrenormal- 
izable local operators. From the high energy viewpoint, there two questions which must be 
answered: 

1. What operators are generated? 

2. With what coefficients? 

For example, the penguin diagrams pictured in Fig. ^ generate, among others, the operators 

C 7 7 = CM (sa^bUF^, 

C 8 8 = C s (fi)(rsYb) L i 1 ,£, (1) 
C 9 9 = CM {s^b) L h^i. 

Perturbative QCD corrections are included by dressing the graphs in Fig. [5] with gluons. The 
leading logarithm approximation, which resums all terms of the form ln n (n/Mw), suf- 

fers from a strong ambiguity in the choice of renormalization scale fi. The resolution of this 
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FIG. 3. Penguin-induced decays of the b quark. 




FIG. 4. Charmless semileptonic b quark decay. 

ambiguity will only properly be resolved by a full next-to-leading order calculation. The 
present state of the art for the coefficient CV(//) is summarized in Ref. ]IJ. The calcula- 
tion is complete to order a s , and partially complete at next-to- leading order. Varying the 
renormalization scale n from m&/2 to mj, one finds a residual scale-dependence uncertainty 
of approximately ±15%. Since O7 is the operator which is primarily responsible for the 
rare decay B — > X s 7, there is a corresponding uncertainty of at least ±30% in the predic- 
tion of this decay rate in the standard model. The coefficients C 8 (//) and C 9 (fi), which are 
responsible for the decay B — > X s £ + £~ , are known with similar accuracy. 

There are also charmless weak b decays, which are rare because their rates are suppressed 
compared to the dominant weak decay mode by the factor \V u b/V c b\ 2 ~ 1CT 2 . Charmless 
semileptonic decays, shown in Fig. |], arise from operators of the form 

76)6^(1 -T5>. (2) 

Since this operator may be written, up to weak and electromagnetic corrections and fermion 
masses, as a product of conserved currents, the coefficient A{jj) suffers from no scale ambi- 
guity. It has been computed to order a s (mb), and the residual uncertainty is small. 

The same is not true of charmless nonleptonic decays, mediated by operators such as 
shown in Fig. |5]. At low energies, these diagrams induce four-quark operators of the form 

Ci Ox = Ci(/i) Mi7 M (l - 7 5 )&i qijl^l - 75)<?2j , 

C 2 2 = C*(ji) Ui^{l - 75)^ <?i,7„(l - 7 5 )<? 2l , (3) 

(4) 
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FIG. 6. New physics contributions to the rare process b — > X s ^/. 

where the indices % and j indicate sums over colors. These operators are not products of 
currents, and they receive renormalizations from perturbative QCD which are as large as 
those received by penguin operators. As summarized in Ref. 0, the coefficients Ci(/x) and 
C2(/i) have now been computed at next-to- leading order. The residual uncertainty, largely 
arising from scheme-dependence, is about ±15%. 

New physics at high energies can also contribute to the coefficients Cj(/i). For example, 
as shown in Fig. || supersymmetric particles, extra scalars, and anomalous trilinear gauge 
couplings can all modify CV(/i) as compared to the standard model. The size of these new 
contributions depends, of course, on the particular model involved. However, because of the 
uncertainties inherent in the perturbative corrections to the standard model, new physics 
will only be observable in rare b decays if causes deviations from the standard model at 
significantly more than the 15% level. This is the most important lesson to be taken from 
the high energy point of view. 



LOW ENERGY VIEWPOINT 

At low energies, we start with an interaction Lagrangian density which is a sum over 
nonrenormalizable operators with coefficients determined at high energies, 

£ = E^)^)' (5) 

The matrix elements of the operators Oi are defined so as to cancel the \i dependence of any 
physical observable. The operators and their coefficients are renormalized at a low energy 
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scale jj, ~ mb, and nonrenormalizable terms are suppressed by powers of the scale M at 
which the interactions become nonlocal and new physics comes into play. 

The challenge, at low energies, is to use the Lagrangian (|5]) to make physical predictions. 
One option is to try to predict exclusive decay modes, such as B — > K*^ or B — > plv. 
However, the theoretical methods available are not entirely satisfactory: the Heavy Quark 
Effective Theory, so useful for b — > c transitions [|3|, is of limited applicability here with 
only light quarks in the final state. Lattice calculations eventually may provide important 
information on exclusive matrix elements, but that is for the most part still in the future. For 
now, one is left to rely on phenomenological models, which for all their occasional successes 
do not provide any controlled approximation to QCD. 

Alternatively, one may consider inclusive decay modes, such as B — > X s 7 or B — > X u tv. 
There has been considerable recent progress [|J in the computation of such quantities in a 
simultaneous expansion in l/m£ and a s {irib) n - It has recently been understood, as well, that 
there are important limitations to such calculations. We shall now review this situation in 
some detail. 

The theoretical analysis of inclusive B decays relies on the Operator Product Expansion 
and perturbative QCD. The partial width V for an operator O to mediate the decay of a B 
to any final state X with the correct quantum numbers is proportional to the square of the 
matrix element, summed over the possible final states, 

r~£|(x|0|i?>| 2 . (6) 

x 

By the Optical Theorem, V may be rewritten as the imaginary part of a forward scattering 
amplitude, 

r ~ Im (B\T{0, O^} \B) , (7) 

which is then expanded simultaneously in powers of a s (jrib) and l/m^. One obtains expres- 
sions for the inclusive partial widths; for example |I]f| . 

T{B - X sl ) oc m\ \C 7 (fi)\ 2 jl + Al ^ A 2 + K(ji) a s (m b ) + ...}, 

T(B - XJLv) cx ml \V ub \ 2 {l + + K*.i. <rn b ) + ...}. (8) 

Here the nonperturbative parameters Ai and A2 are defined by hadronic matrix elements [Q, 

Ai = (B\ b(iD) 2 b \B)/2m b , 

A 2 = (B\ b(-i a ^)G^b \B)/2m b . (9) 

It is straightforward to find K s j 4 = — vr 2 ) 0, while the perturbative correction K{fi) 

is too messy to be illuminating [|TJ. 

There are a number of sources of uncertainty in the expressions The operator 

b(—^a' JlU )G^ u b violates the Heavy Quark Spin Symmetry and may be measured directly 
from the B-B* mass difference, A2 ~ 0.12 GeV 2 . However, Ai can not be measured directly; 
instead, one must rely on phenomenological models. While this is unfortunate, if we assume 
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FIG. 7. A typical one-loop radiative correction to T{0,0^}. 
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FIG. 8. The gluon propagator replaced by the sum of self-energy graphs. 



Ai < 1 GeV 2 , then inspection of Eq. (§) shows that the error induced in the partial widths 
T is 10% or less. Higher order nonperturbative corrections, of order l/m\, are expected to 
be at the level of a few percent. 

The primary sources of uncertainty in Eq. (|8]) are the value to take for the bottom mass 
rrib, and uncomputed higher order radiative corrections. Because of the overall factor of 
m\, the theoretical partial widths are extremely sensitive to this parameter. For example, 
allowing m^ to vary over the range 4.5 GeV < m^ < 5.0 GeV induces an uncertainty in V of 
approximately 50%. While lower values for m b seem currently to be preferred, the issue is 
still quite unsettled.^] 



Higher order radiative corrections to the partial widths have recently been considered 
by a number of authors. In particular, attention has been paid to a set of corrections 
which are dominant in the limit of large Nf (large number of quark flavours), and which 
in the real world still may be particularly large. These come from taking the one loop 
radiative correction to the time-ordered product (0), an example of which is shown in Fig. [7], 
and replacing the gluon propagator with a sum of self-energy bubbles. If this replacement, 
which is illustrated in Fig. || is carried out to all orders and then extrapolated to the physical 
Nf, it amounts to replacing the strong coupling constant a s {m b ) by its running value a s (p 2 ) 
evaluated at the loop momentum. 

The BLM scale-setting prescription |10[ requires that one perform the substitution shown 
in Fig. [8] to leading order; the two-loop contribution to the radiative correction is then 
expected on general grounds to be parametrically large. Once this part of the two-loop 



There is an ongoing controversy over issues as fundamental as the proper definition of mj. I will 
not review this discussion here. 



A. Higher order radiative corrections 
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FIG. 9. A cartoon of the lepton energy spectrum in B — > X u lv, along with the kinematic cut. 

computation has been done, one adjusts the scale \i in the one-loop result to absorb it. A 
recent application of this criterion to the inclusive rate for B — > X u lv indicates that the 
appropriate scale for this process is \i ~ mb/ 10 rather than nib [TT|. A more complicated 
scale-setting procedure which resums all orders in the bubble sum (but which suffers from 
a certain lack of uniqueness) does not, in general, lead to quite such a low value of fi ]P2| |. 
although the two-loop corrections are, of course, still quite large. 

This treatment of higher order radiative corrections leaves us with two questions. 

1. Should such a low renormalization scale be taken seriously? If so, then clearly the entire 
program of computing inclusive rates perturbatively is in trouble. If not, then one still has 
to do deal with the fact that two-loop corrections are much larger than one might naively 
have thought. 

2. If there is a class of diagrams which is unusually large, can perturbation theory be 
improved in a sensible way? Once such a resummation has been performed, can one show 
that the remaining uncertainties are likely to be small? 

B. The need for endpoint spectra 

Final states with charm present an enormous background to rare B decays. For example, 
the decay B — ► XJLv obscures B — > X u £v, and B — > Dtt° — ► D77 presents a problematic 
background to B — > X s ^. Typically, strict kinematic cuts are used to exclude such process. 
For example, studies of rare decays accept only leptons and photons with energies in the 
range 2.2 GeV < Ej»,E 1 < 2.7 GeV, beyond the kinematic endpoint for charm in the final 
state. Hence it is necessary for theorists to compute not only partial widths T, but inclusive 
lepton spectra dT/dE within 20% or so of the endpoint. A cartoon of a lepton energy 
spectrum, along with the kinematic cut, is shown in Fig. |9|. 
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The theoretical problem is that the OPE does not converge when resticted to the lepton 
or photon energy endpoint. What is computable is not the full differential spectrum dT/dE 
but rather moments of this spectrum, obtained by weighting the differential spectrum by 
some function and then integrating. Introducing the scaled energy variable y = 2Eg^/mf 1 , we 
thus "smear" with a weighting function with support only in the small region 1 — 5 < y < 1. 
The size 5 of the smearing region controls the convergence of the OPE. For 5 ~ Aqcd /^rib ~ 
10%, all orders in the 1/rrib expansion contribute equally |J, and the leading terms in the 
expansion (||) is clearly insufficient. Of course, this is only an order of magnitude estimate, 
and how the OPE converges for the experimentally chosen upper value of 5 cannot be 
determined from such general considerations. At this point, then, a certain amount of faith 
is required in the interpretation of the smeared theoretical spectra. 

An interesting by-product of this analysis is the result that the same infinite sum of 
terms in the expansion determines the shape of the endpoint spectrum in B — > X s 7 

and in B — > X u tv H0]. Whether this relation yields useful predictive power is still to be 
seen. 



C. Sudakov Logarithms 



Another source of uncertainty in the shape of the endpoint spectrum comes from Sudakov 
logarithms |13| . For example, the perturbative corrections to the lepton energy spectrum in 
B — > X u tv is extremely singular near the endpoint y = 1 fl4|j: 



dr 

dy 



dTo 
dy 



2as 
3vr 



31 

\n 2 (l-y) + —\n(l-y) 



(10) 



where dT /dy is the spectrum at tree level. At order a 2 s) the leading singularity is ln 4 (l — y), 
and so forth. These Sudakov double logarithms may be resummed into an exponential 
suppression factor: 



dr _ dr 

dy dy 6XP 



2a, 
3vr 



ln 2 (l-y) +... • 



(11) 



This leading behaviour is actually stronger very near y — 1 than that given by the nonper- 
turbative power corrections, but it is calculable. 

What must be suppressed are the leading uncalculated corrections, which is accomplished 
by smearing over a region 5 large enough that they may be neglected. In the large mb limit, 
this requires that we smear over a region formally much larger than 8 ~ Aqcd/^6, given by 
the condition 



S > exp | - '^njajjnb) j 



(12) 



In this strict limit, then, all nonperturbative corrections to the endpoint shape would be 
irrelevant. But for realistic w 4.8 GeV, and the given experimental smearing region 
5 ~ 0.1 ~ 0.2, do the uncalculated Sudakov effects actually dominate the nonperturbative 
power corrections? It is difficult to guess, based only on naive power counting arguments. 
Explicit calculations of the subleading Sudakov logarithms may help clarify the situation. 
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FIG. 10. The one instanton contribution to T{0,0^}. 

For technical reasons, the Sudakov corrections to the smeared photon spectrum in B — ► 
Xg'y are under much better theoretical control than in B — > X u £u || . They do not introduce 
unmanageable uncertainties into the computation of the weighted spectra. 



D. Instantons 

Finally, there are possibly large contributions to energy endpoint spectra from instantons. 
These arise because the light quark which is produced in the short-distance interactions can 
propagate in an instanton background, as pictured in Fig. [II]. Chay and Rey computed, 
in the dilute instanton gas approximation, the one instanton contribution to dT/dy for 
B — > X u lv and B — > [|nj. Their result diverges dramatically at the endpoint, as y — > 1. 



The contribution to B — > X s ^ is nonetheless small and under control when one computes 
weighted spectra, but the same is not true for B — > X u iv. Instead, one finds that the one 
instanton contribution is entirely untrustworthy in the experimentally defined window. 

The one instanton contribution goes bad in this region presumably because multi- 
instanton configurations begin to be important. We have used the one instanton calculation 



as the motivation for a crude ansatz for the multi-instanton result in this region [|16j. This 
ansatz incorporates, as much as possible, the reliable information from the one-instanton 
calculation. When we vary this naive "best guess" ansatz by two orders of magnitude, we 
find that over most of the ansatz parameter space, the instantons do in fact dominate the 
weighted endpoint spectra. 

One must be careful about interpreting this result. It is potentially interesting only in 
a negative sense. On the one hand, the actual numbers certainly cannot be believed; by 
no means do we claim to have computed the correct multi-instanton contribution. On the 
other, we have failed to find any justification for ignoring the instantons in the endpoint 
region. In light of this equivocal situation, one may well wonder whether one can still trust 
the relationship between the endpoint spectra for B — > X s ^ and B — > X u iv proposed in 
Refs. ||[7|. This is a situation badly in need of clarification. Invocations of faith, one way 
or the other, will not be sufficient; a more sophisticated estimate of instanton contributions 
is what is required. Such an estimate could show, for example, that our ansatz for the 
multi-instanton contribution is entirely too crude, and that other techniques can be used to 
prove that the multi-instanton contribution is necessarily negligible. We certainly hope that 
this will prove to be the case. 
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CONCLUSIONS 



We may summarize the status of the theory of rare B decays from each of our two 
viewpoints: 

Low Energy Viewpoint: 

1. Exclusive decay rates are extremely difficult to compute reliably. One must resort to 
models and other uncontrolled assumptions, a situation which is most unsatisfactory. 

2. Inclusive calculations, by contrast, may be performed in a controlled expansion in powers 
of a s and l/ra&. However, there remain unresolved uncertainties about 

a. uncomputed higher orders in a s and the renormalization scale //; 

b. Sudakov double logarithms near the lepton energy endpoint; 

c. instanton contributions near the lepton energy endpoint. 

3. The decay B — > X s 7 is in much better shape with respect to Sudakov and instanton 
corrections than is B — > X u £u. Hence, while the calculation of B — > X s 7 is itself perhaps 
fairly secure, the proposed relationship between the endpoint spectra in B — > X s ^ and 
B — > X u £v may well be threatened by these effects. 

High Energy Viewpoint: 

1. In view of the significant uncertainties in existing theoretical calculations, only modifica- 
tions to the Standard Model which affect rare decays at the 50% level or higher are likely to 
be experimentally detectable. Small modifications, say at the 10% level, are unlikely ever 
to be seen. 

2. There is considerable room for the situation to improve, and much work remains to be 
done. 
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